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5WFr OBSERVATIONS OF THE 2006 OUTBURST OF THE RECURRENT NOVA RS OPHIUCHI: 

III. X-RAY SPECTRAL MODELLING 

N. M. H. Vaytet'-^, t. J. 0'BRIEN^ K. L. Page"*, M. F. Bode^, M. Lloyd'' and A. P. Beardmore^ 

ABSTRACT 

Following the Swift X-ray observations of the 2006 outburst of the recurrent nova RS Ophiuchi, we developed 
hydrodynamical models of mass ejection from which the forward shock velocities were used to estimate the 
ejecta mass and velocity. In order to further constrain our model parameters, here we present synthetic X- 
ray spectra from our hydrodynamical calculations which we compare to the Swift data. An extensive set of 
simulations was carried out to find a model which best fits the spectra up to 100 days after outburst. We find a 
good fit at high energies but require additional absorption to match the low energy emission. We estimate the 
ejecta mass to be in the range (2 — 5) x 10^^ Mg, and the ejection velocity to be greater than 6000 kms^' (and 
probably closer to 10,000 kms^'). We also find that estimates of shock velocity derived from gas temperatures 
via standard model fits to the X-ray spectra are much lower than the true shock velocities. 

Subject headings: hydrodynamics — shock waves — novae, cataclysmic variables — stars: winds, outflows 
— stars: individual (RS Ophiuchi) — X-rays: binaries 



1. INTRODUCTION 

RS Ophiuchi (RS Oph) is one of th e mos t well-studied of 
the ten known recurrent novae (RNe; lAnupa ma 2008). The 
central system is a b inary comprising a white d warf (WD) 
and a red giant (RG) dDobrzvcka & Kenvonll 19941) . A classi- 
cal nova outburst involves mass transfer from the companion 
to the WD. The build-up of pressure and temperature in the 
degenerate layer of accreted hydrogen eventually leads to a 
thermonuclear runaway (TNR) on the surface of the WD, re- 
sulting in the high-speed eje ction of a shell of material into 
the circumstellar medium (St arrfieldl 120081) . Novae which 
have been observed in outburst more than once are called 
RNe. The shock interaction of the ejecta with the surround- 
ing medium has been found to heat the gas to temperatures of 
. 10^ - 10^ K, yielding hard X-rav radiation (iLloyd et al.l lT992l: 
■ lO'Brien. Llovd & Bodd[T994l: iMukai & Ishidall2001[) . In the 
\ case of RS Oph, the ejecta run into the surrounding dense RG 
wind. Soft X-ray emission is also expected to be revealed later 
in the outburst from a centra l WD close to Eddington lumi- 
\ nosity (iKrautter et al.l Il996t iBalman. Krautter & Oegelmaiil 
. 119981) as is seen in RS Oph (see below). 

RS Oph has h ad outbursts re c orded in 1898, 1933, 1958, 
' 1967, 1985 (see Rosino' 119871 |Rosino & liiima 1987|) and 
most recently on 2006 February 12 (Narumi et al. 20 06i), with 
' possible additional outbursts in 1907 and 1945 dSchaefed 
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120041 lOppenh eimer & Mattei|[T993h . Its binary system ha s 
an orbital period of 455.72 ±0.83 days dFekel et al.ir2000h . 
The mass of the WD is suspected to be close to the Chan- 
drasekhar limit from various considerations, including the 
speed class of outburst an d the relatively shor t recurrence time 
between outbursts (e.g. ' Yaron et aTl I2005L ISokoloski et alj 
2006, Hachisu et al. 200 7). The distance to RS Oph is 1.6 
kpc (Bode 1987; see also Barrv et al.ll2008h . 

RS Oph was detected in X-rays just 2 days after the 2006 
outburst with the RXTE satellite (ISokoloski et al.ll2006 *) and 
3.17 days after ou tburst using the XRT instrument on board 
the Swift satellite dBode et al.ll2"006l hereafter |P^E1- The 
Swift observations were performed on average every one or 
two days during the first 100 days, covering all the phases 
of the remnant's evolution. The rate of observations sub- 
sequently decreased but continued for more than a year 
(Page et al. 2008). RS Oph was seen initially as a source of 
fairly hard X-rays, gradually softening with time. Around 26 
days after outburst, the very bright and soft component from 
the Super-Soft Source (SSS) phase of nuclear burning on the 
surface of theJWD appeared in the spectrum (Osborne et al] 
2006 1 : iHachisu et all BooH iPage et al.l l2008i: lOsborne etid] 

In an effort to improve on the models of the 1985 outburst 
and, in particular, to address the new observations of the early 
phases of its evoluti on, we developed new hydro dynamical 
models for RS Oph (IVaytet. O'Brien & Bode"2007. hereafter 
Paper II). We u sed a one-dimension al Eulerian second order 
Godunov code dO'Brien et al. 1994) to simulate the outburst 
where mass-loss in the form of a fast wind from the WD ran 
into a surrounding slow RG wind. The scheme took into ac- 
count all three phases of the remnant's evolution, and allowed 
us, after a thorough analysis of the shock velocities, to de- 
termine that the ejecta had a low mass ('-^ 2 x 10^^ Mg) and 
a high ejection velocity. The models also showed that radia- 
tive cooling had an important impact on the evolution of the 
shocks; the amount of energy radiated by the hot ejecta was 
correlated to the rate at which the forward shock velocity de- 
cre a sed in tim e . 

In iPaper III our analysis reli ed on sin gle-temperature fits to 
the Swift X-ray spectra (from iPaper J) in order to determine 
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post-shock temperatures and hence the shock velocities. Here 
we improve on this approach by calculating full X-ray spec- 
tra directly from our hydrodynamical results using XSPEC0 
and by including the effects of overlying absorption and in- 
strument response. Results from an extensive set of simula- 
tions are compared to the Swift observations for the first 100 
days after outburst in order to obtain best-fitting values for the 
model parameters. 

2. 5W/FrXRT OBSERVATIONS 

Some 70 observational epochs of RS Ophiuchi were ob- 
tained using the X-ray Telescope (XRT) onboard the Swift 
observatory during the first 100 days after the 2006 outburst. 
As in Paper I, the data were processed using the standard 
Swift software and source spectra were extracted from the 
cleaned Windowed Timing (WT) mode (f < 90 days) and Pho- 
ton Counting (PC) mode (f > 90 days) event lists. Because 
of the high count rate of the X-ray source, most of the WT 
were piled-up. Thus, the core of the PSF was excluded: an 
outer box of width 60 x 20 pixels was used, with the cen- 
tral 10 pixels (1 pixel = 2"36) removed. Similarly, the PC 
data were piled-up: here, an annulus with outer radius 30 
pixels was used, with the inner exclusion radius decreasing 
from six to zero pixels between days 90 and 100. Background 
spectra were obtained from nearby source-free regions, us- 
ing a 60 X 20 pixel box for the WT data and a 60 pixel ra- 
dius circle for the PC. We corrected for the fractional loss 
of the point-spread function caused both by the exclusion of 
the central (piled-up) region and the positioning of the source 
over the bad CCD colur nns (caused by mic rometeorite dam- 
age on 2005 May 2007: lAbbev et al.ll2006l) when calculating 
the count rate and flux. 

Grades 0-2 and 0-12 were chosen for the spectral analysis 
of the WT and PC data respectively. The FTOOL xrtmkarf 
was used to generate suitable ancillary response function files, 
and these were used in conjunction with the relevant response 
matrices. All spectra were grouped to a minimum of 20 counts 
per bin in order to facilitate the use of statistics in XSPEC, 
the results of which are shown in Fig.[T](black crosses) for 8 
epochs. 

3. CALCULATION OF SYNTHETIC X-RAY SPECTRA 

Our hydrodynamical models from lPaper Ill predict the phys- 
ical conditions of the hot shocked gas at any time during the 
outburst. In practise the model comprises many concentric 
spherical shells each with a different temperature and density. 
Since it would be too computationally intensive to calculate 
the X-ray emission from each shell separately, we grouped 
them into 30 logarithmically-divided temperature bins from 
10"* to 10*^ K, weighting each shell's contribution by its emis- 
sion measure (experiment suggested 30 bins were sufficient 
to give an accurate result). These were then fed as a 30- 
component plasma into XSPEC (version 12), an X-ray data re- 
duction and spectrum calculation code. This was used to com- 
pute the X-ray spectrum with the MEKAL model for emission 
from hot, diffuse gas, including line emissions from C, N, O, 
Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe and Ni, along with free- 
fre e, free-b ound and two-photon e mission (for more details 
see |\^ytet 2009'; iMewe et al.| [T98g). 

The interstellar equivale nt hydrogen absorption column for 
our models was taken from lHiellming et aTl (119861) who found 

^ http://swift.gsfc.nasa.gov/docs/software/lheasoft 



nn = (2.4±0.6) x 10^' cm^^. We have also included a contri- 
bution from the circumstellar absorption due to the RG wind 
surrounding the hot, shocked gas behind the forward shock. 
At each epoch we calculate an average value for this circum- 
stellar absorption across the face of the sphere of shocked gas. 

Of course, only the unshocked neutral part of the RG wind 
will contribute to absorption. It is however non-trivial to know 
the ionisation state of the RG wind at all times, since a large 
portion of it could be ionised by the radiation emitted dur- 
ing the outburst. This requires not only knowledge of the re- 
combination rates in the RG wind but also knowledge of the 
luminosity of the central radiation source at all times, which 
is clearly not const a nt nor simple to model (see f or exam- 
ple lPage et all2008h . lO'Brien. Bode & KahnI ( 11992b assumed 
that the whole RG wind was flash-ionised by UV radiation at 
the very beginning of the outburst, then computed recombina- 
tion rates inside the RG wind. Here we simply use upper and 
lower extremes; either the RG wind is fully neutral or fully 
ionised. We are then confident that the correct absorption has 
to lie between these limits since they do not depend on any 
assumptions about the central radiation field. 

The final step in the creation of the synthetic spectra was to 
fold the results through the Swift XRT instrumental response 
matrix. The response used here was the same as for the ob- 
servations mentioned above, with the corresponding ancillary 
file specific to each observation to provide the precise effec- 
tive area of the detector and correct for bad channels. We used 
the canonical distance of 1.6 kpc for normalisation. 

4. RESULTS 

4.1. Initial simulations from Waper lA 

Mode l runs 1 to 16 presented in this section are identical to 
those in iPaper if . see Table [1] for a list of parameter values. 
Run 14 is not included since it represented a continuous fast 
wind ejection during the entire duration of the simulation, thus 
rendering parameters such as ejected mass Mgj and outburst 
energy Eq meaningless. 

Fig. [T] shows the synthetic spectra for five of these sim- 
ulations compared with the Swift observations. Results are 
shown in separate panels for days 3, 5, 11, 18, 29, 50, 75 
and 100 after outburst. Each panel includes results assuming 
the overlying RG wind is either fully neutral (solid line) or 
fully ionised (dashed line). We first consider the results from 
run 2 (shown in black), a fairly typical example. At early 
times (3 — 18 days), before the SSS component from the WD 
remnant appears at the soft end of the spectrum, the synthetic 
spectra clearly do not match the observations; the model is 
overestimating the soft X-ray count rates and underestimating 
the hard count rates. We do not include the SSS component 
in our models so for later times (day 29 and after), the Swift 
data at the soft end of the spectrum has been omitted (energies 
< 1 keV). At these late times, at the hard end of the spectrum 
the model count rates are very closely matched to the obser- 
vations, suggesting that the hard count rates have decreased 
faster in the observations than in the model. 

In order to link the synthetic spectra with the hydrodynam- 
ics, we recall Fig. 5 in Paper II where we concluded that the 
deceleration rate of the forward shock was an important as- 
pect of the shock models. We concluded that runs 3 and 1 1 
have forward shock deceleration rates lower and higher than 
run 2, respectively. Their X-ray spectra are shown in Fig. [T] 
(orange for run 3 and blue for run 11). 

In run 3 (orange), the increased Mej increases the X-ray 



X-RAY SPECTRAL MODELLING OF RS OPH (2006) 



3 



o 0.5 1 2 5 10 0.5 1 2 5 10 




0.5 1 2 5 10 0.5 1 2 5 10 

Energy (keV) Energy (keV) 



Figure 1. Synthetic X-ray spectra for runs 2 (black), 3 (orange), 5 (green), 1 1 (blue) and 13 (dark pink) at days 3, 5, 11, 18, 29, 50, 75 and 100. The solid and 
dashed lines represent a fully neutral and fully ionised RG wind, respectively. The black crosses represent the Swift data. The Swift data at low energies for days 
29 and after showing the SSS component have been removed since no SSS is included in our models. 



count rates, this effect being strongest at the hard end of the 
spectrum. The hard X-ray count rates decrease more slowly 
than in run 2, which is due to the higher inertia that the shell 
carries due to its higher mass. This is best seen at days 18 
to 50. Days 75 and 100 should be ignored in this case as by 
this time the shoc k has reached the edge of the RG wind; see 
Fig. 4 in lPaper 11 

In the case of run II (blue), the high fast wind velocity 



V2 also increases the X-ray count rates, but this time over 
the whole energy spectrum, most probably simply due to the 
higher outburst energy compared to run 2. In run 1 1 both hard 
and soft X-ray count rates decrease much faster than in runs 2 
and 3, falling to the bottom of the visible plot area by day 75. 
Note that for run 1 1 the shock does not reach the edge of the 
wind on a 100-day timescale. Finally, the impact of absorp- 
tion in the overlying high-density neutral RG wind in run 1 1 
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(solid blue line) is clearly visible on the soft part of the spec- 
trum at early times where the predicted count rates are more 
consistent with observations. 

Further, by simply increasing V2 and leaving Mgj untouched 
as in run 5 (Fig-B green lines), we note that the count rates are 
also greatly increased across the spectrum. Runs 3 and 5 have 
identical outburst energies and RG wind densities; it is hence 
not surprising to see that their X-ray spectra are very similar 
at late times (day 29 and after) apart from the hard count rates 
being slightly higher in run 3 at day 50. Comparing run 5 
with run 1 1 reveals that increasing V2 can only be effective in 
enhancing hard count rates at early times if Mej is not too low. 

We conclude that, as one might expect, greater decelera- 
tion of the forward shock results in faster reduction in hard 
X-ray count rates. This is also sensitive to the ejected mass 
Mej; a high-mass shell will carry more momentum and there- 
fore decelerate more slowly. A high ejection velocity (V2 > 
5000 km s^^) greatly overestimates the soft count rates. Note 
that the strong Fe XXV K line at 6.7 keV appears wider in 
the Swift observations than in our models (run 2 for example). 
We do not include Doppler broadening in our models but it 
should also be noted that a second Fe XXVI line (6.97 keV) 
appears very close to the first in some of our models (runs 3, 



5 and 13) which would not be resolved by the XRT and would 
broaden the appearance of the Fe line. Finally, it is clear from 
runs 2 and 3 that a relatively low-density (even entirely neu- 
tral) RG wind appears not to produce sufficient absorption of 
low- energy X -rays at early times. 

In iPap er II run 13 provided the best fit to the forward 
shock velocities. However the extremely high implied V2 of 
16,300 km s^', for which there was little evidence elsewhere, 
cast some doubt on the validity of the model. The synthetic 
spectra for run 13 are shown in Fig. [T](dark pink). These 
clearly show run 13 does not match the observations of the 
2006 outburst of RS Oph. Although the desired large amount 
of absorption is present at early times, the count rates are over- 
estimated across the spectrum at every epoch. 

4.2. Fitting synthetic spectra to the Swift data 

A fuller exploration of the parameter space of the hydro- 
dynamical model has now been undertaken by carrying out 
a new extensive set of simulations using the same code as in 
Paper II. This brings the total number of runs to 88. All the 
parameters for runs 1 to 88 are listed in Table [T] (run 14 ex- 
cluded for reasons given above). 



Table 1 ASPHERE CODE PARAMETER SPACE 



Pun 


RG wind 


Fast wind 


Fast wind 


Outburst 


Ejected 


Ko wmu 


IVIass ratio 


Ko winu 


Fast wind 


Mean flux difference 




nitiss-loss 


nitiss-loss 


velocity 


energy 


mass 


mass 






phase 


total 


early 


late 




M\ 


M2 


V2 


Eo 


Mej 


MRGW 


R 


M\ 1 It 


'o 




<J>£ 




(> 


;10"' M(.-,/yr) (> 


;10"^ M .,/yr) 


(kms"') 


{xW^^ erg) 


(X 10"'' Mq) 


(xlO^* M,>} 




(xlO'2 g/cm) 


(days) 




(f<30d) (f>30d) 


1 


1.43 


158.59 


3000 


1.85 


0.24 


3.00 


0.080 


6.00 


7.0 


0.961 


0.940 


0.965 


2 


1.43 


730.62 


3000 


8.62 


1.10 


3.00 


0.367 


6.00 


7.0 


0.574 


0.748 


0.537 


3 


1.43 


3171.80 


3000 


36.94 


4.78 


3.00 


1.595 


6.00 


7.0 


0.484 


0.570 


0.465 


4 


1.43 


730.62 


1398 


1.85 


1.10 


3.00 


0.367 


6.00 


7.0 


0.997 


0.995 


0.998 


5 


1.43 


730.62 


6251 


36.94 


1.10 


3.00 


0.367 


6.00 


7.0 


0.611 


0.311 


0.676 


6 


1.43 


730.62 


3500 


11.58 


1.10 


3.00 


0.367 


6.00 


7.0 


0.442 


0.661 


0.395 


7 


1.43 


730.62 


4000 


15.13 


1.10 


3.00 


0.367 


6.00 


7.0 


0.411 


0.567 


0.378 


8 


1.43 


5113.73 


3000 


8.62 


1.10 


3.00 


0.367 


6.00 


1.0 


0.528 


0.628 


0.506 


9 


1.43 


340.97 


3000 


8.62 


1.10 


3.00 


0.367 


6.00 


15.0 


0.660 


0.746 


0.642 


10 


6.34 


730.62 


3000 


8.62 


1.10 


13.32 


0.083 


26.67 


7.0 


0.995 


0.974 


1.000 


11 


6.34 


79.30 


12000 


14.77 


0.12 


13.32 


0.009 


26.67 


7.0 


0.853 


0.259 


0.980 


12 


1.43 


1148.03 


3969 


8.62 


1.10 


3.00 


0.367 


6.00 


7.0 


0.567 


0.697 


0.539 


13 


9.52 


158.59 


16300 


46.73 


0.20 


19.98 


0.010 


40.00 


6.0 


17.956 


6.491 


20.413 


15 


1.43 


298.62 


3000 


3.48 


0.45 


3.00 


0.150 


6.00 


7.0 


0.634 


0.849 


0.588 


16 


1.43 


497.66 


3000 


5.80 


0.75 


3.00 


0.250 


6.00 


7.0 


0.620 


0.786 


0.584 


17 


7.93 


158.59 


9000 


9.50 


0.14 


16.65 


0.008 


33.33 


4.0 


0.986 


0.921 


1.000 


18 


52.33 


852.29 


4500 


19.14 


1.10 


110.00 


0.010 


220.20 


6.0 


4.160 


0.728 


4.895 


19 


7.93 


774.82 


8000 


55.00 


1.00 


16.65 


0.060 


33.33 


6.0 


19.847 


6.348 


22.740 


20 


6.34 


126.87 


10000 


16.42 


0.19 


13.32 


0.014 


26.67 


7.0 


0.748 


0.375 


0.828 


21 


6.34 


174.45 


7000 


11.06 


0.26 


13.32 


0.020 


26.67 


7.0 


0.975 


0.863 


1.000 


22 


7.93 


451.57 


7000 


24.54 


0.58 


16.65 


0.035 


33.33 


6.0 


1.099 


1.276 


1.061 


23 


7.93 


323.49 


8000 


22.96 


0.42 


16.65 


0.025 


33.33 


6.0 


0.908 


0.877 


0.915 


24 


7.93 


646.97 


7000 


35.16 


0.83 


16.65 


0.050 


33.33 


6.0 


10.241 


3.533 


11.678 


25 


7.93 


646.97 


6000 


25.83 


0.83 


16.65 


0.050 


33.33 


6.0 


1.223 


1.534 


1.157 


26 


7.93 


646.97 


5000 


17.94 


0.83 


16.65 


0.050 


33.33 


6.0 


0.939 


0.655 


0.999 


27 


7.93 


774.82 


4500 


17.40 


1.00 


16.65 


0.060 


33.33 


6.0 


0.958 


0.763 


0.999 


28 


7.93 


1035.15 


4500 


23.25 


1.34 


16.65 


0.080 


33.33 


6.0 


0.775 


0.408 


0.853 


29 


7.93 


1293.95 


4500 


29.06 


1.67 


16.65 


0.100 


33.33 


6.0 


2.973 


1.578 


3.272 


30 


7.93 


1164.56 


4000 


20.67 


1.50 


16.65 


0.090 


33.33 


6.0 


0.901 


0.466 


0.995 


31 


7.93 


1578.30 


4500 


35.45 


2.04 


16.65 


0.122 


33.33 


6.0 


8.992 


2.471 


10.390 


32 


7.93 


1940.92 


4500 


43.59 


2.51 


16.65 


0.150 


33.33 


6.0 


14.352 


3.447 


16.688 


33 


7.93 


2587.89 


4500 


58.12 


3.34 


16.65 


0.201 


33.33 


6.0 


17.011 


4.900 


19.606 


34 


7.93 


1578.30 


4000 


28.01 


2.04 


16.65 


0.122 


33.33 


6.0 


2.124 


1.351 


2.289 


35 


7.93 


1940.92 


4000 


34.44 


2.51 


16.65 


0.150 


33.33 


6.0 


8.035 


2.174 


9.291 


36 


7.93 


2587.89 


4000 


45.92 


3.34 


16.65 


0.201 


33.33 


6.0 


14.954 


3.296 


17.452 


37 


7.93 


1578.30 


3500 


21.44 


2.04 


16.65 


0.122 


33.33 


6.0 


0.862 


0.368 


0.968 


38 


7.93 


1940.92 


3500 


26.37 


2.51 


16.65 


0.150 


33.33 


6.0 


0.986 


0.915 


1.001 


39 


7.93 


2587.89 


3500 


35.16 


3.34 


16.65 


0.201 


33.33 


6.0 


7.962 


2.065 


9.225 


40 


7.93 


155.28 


12000 


24.80 


0.20 


16.65 


0.012 


33.33 


6.0 


1.366 


1.745 


1.284 


41 


7.93 


182.39 


12000 


29.13 


0.24 


16.65 


0.014 


33.33 


6.0 


4.023 


2.979 


4.247 


42 


7.93 


232.91 


12000 


37.20 


0.30 


16.65 


0.018 


33.33 


6.0 


13.453 


4.468 


15.378 


43 


7.93 


284.67 


12000 


45.46 


0.37 


16.65 


0.022 


33.33 


6.0 


19.062 


5.674 


21.931 


44 


7.93 


155.28 


10500 


18.99 


0.20 


16.65 


0.012 


33.33 


6.0 


0.885 


0.349 


0.999 
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TABLE 1 - Continued 



Run 


M\ 


Ml 


V2 




M 


M 




M\ ju 


'0 








45 




182 39 


10500 


99 10 
ZZ. JU 


24 


IfS ^^S 

lO.OJ 


014 


11 11 
JJ. J J 




o.u 


S4S 


750 


866 


46 


/ .yj 


9^9 Ql 


10500 


28 48 


10 

U. JU 


16 65 


01 s 

U.Ul O 


11 11 
J J. J J 


fy 
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Notes. -The slow wind velocity is « = 15 kms throughout (PagerlT). Run 14 is omitted (see section l4~l1 and PagerlT). 



In 'Paper IJ, we concluded that the deceleration rate of the 
forward shock was directly related to the amount of energy 
radiated away by the ejecta. This in turn appeared to be related 
to the ratio of the ejected mass to the mass contained in the 
RG wind; the lower the ratio, the higher the fraction of energy 
which was radiated away. In Table. [T] we also list this ratio 
R = Mgj / Mrgw for each run. 

In order to assess the goodness of fit of each run's synthetic 
spectra to the Swift data, we extracted fluxes from our mod- 
els in order to compare them directly to the fluxes observed 
in some 70 Swift observations between day 3 and day 100. 
In all cases, the fluxes have been taken in the energy range 
1 — 10 keV in order to avoid contamination from the SSS 
emission. The absolute value of the difference between the 
model and observed fluxes as a fraction of the observed flux 
was calculated for each epoch. These were then summed over 
all epochs and an average was taken, giving the total mean 
fractional flux difference <Pt- We also define two other quan- 
tites <Pe and <Pl which are the early (f < 30 days) and late 
(30 < f < 100 days) mean fractional flux differences respec- 
tively. These values are listed for all runs in Table [T] All the 
values quoted are for runs with a fully-neutral RG wind. It 



is evident from the Swift observations that some absorption is 
present in the spectra, and listing the mean flux differences 
for the synthetic spectra with an ionised RG wind is therefore 
unecessary. Moreover, the effect of circumstellar absorption 
is highest at the soft end of the energy spectrum (see Fig. [U 
and since we only consider energies in the range 1 — 10 keV, 
total mean flux differences for ionised and neutral RG winds 
are in almost all cases very similar. 

The total mean flux differences are a useful way of identify- 
ing the regions in the parameter space which yield best fitting 
runs, however it is not an absolute measure of the quality of 
the fits. For instance run 1 1 would be classed as one of the best 
fitting runs according to (S>e (Table [U, but the model spectra 
at days 3 and 5 (see Fig.IIJ are clearly not good fits to the data 
as the spectral slope is very different to the observations. The 
over-estimate in medium energy counts (1 .0 — 2.5 keV) com- 
pensates for the underestimate at the hard energy end, giving 
a total flux very close to the observed flux. We exclude any 
runs with large <t> values. We also examined all the spectra 
by eye to ensure that we had not excluded runs which had the 
correct spectral shape but were lying either above or below 
the Swift data at all epochs. Indeed, such a systematic error 
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could be corrected simply by varying the distance to RS Oph , 
which has been subject to some debate dBarrv et alj 120081) . 
After this exhaustive analysis, we have come to the following 
conclusions (ejecta masses are given in 

• For all runs, in order to reproduce the correct amount of 
hard X-ray emission, we require V2 > 6000 kms^' 

• If Mej > 1.5 (x 10^^ Mq), spectra can be a good match 
at early times (f < 5 days) but are much too bright at 
later times (e.g. run 33) 

• If 0.8 < Mej < 1.5: 

- if Ml is high, the X-ray flux is either too high 

(run 19) or fades too quickly (run 25) 

- if Ml is low, a high V2 > 6000 kms^' forces the 

fluxes to be too high and the deceleration rate 
of the forward shock to be too low (run 88) 

• If 0.2 < Mej <0.8: 

- if Ml is high, the X-ray flux is either too low for 

very low masses (run 48) or fades too quickly 
(runs 23, 45) 

- if Ml is low, hard X-ray counts are too low for 

V2 < 8000 km s"^ (run 73) 

- a good fit is found for the hard X-rays at all 

epochs if Ml is low and V2 > 8000 kms^' 

• If Mej < 0.2, the hard X-ray flux is too low and the 

fluxes at all energies fade very rapidly (runs 
11, 17) 

This step by step elimination process leads us to believe that 
the mass ejected in the outburst Mgj is in the range (2 — 8) x 

10^^ Mq. We find our best fitting model to be run 81 (see 
Fig. 121), although even this does show a small overestimate in 
hard X-rays around day 29. Its main characteristics are its low 
Mej = 2.7 X 10-^ Mq and high V2 = 10'^ kms^^ 

The majority of the runs produce too many low energy X- 
rays. We have investigated whether this might be due to insuf- 
ficient circumstellar absorption. In order to naturally increase 
the absorption, we simply increased the slow wind mass-loss 
rate (see for example the difference in the amount of absorp- 
tion in runs 2 and 1 1) but of course this has a significant effect 
on the shock dynamics. With a denser RG wind, if we wish 
to produce shock velocities high enough to be consistent with 
observations, we have to vastly increase V2. This leads to an 
overestimate of the soft and medium energy count rates (see 
run 13 for example). In order to obtain flat spectra at early 
times (days 3 and 5, say), we have had to keep V2 in the range 
^ 5000— 10000 kms^' and consequently the mass-loss rate 
into the RG wind Mi also needs to be relatively low. A mid- 
range Mej was then required to obtain the correct hard X-ray 
decay rate, as in run 8 1 . 

An alternative way of addressing the deficit in absorption 
of soft X-rays is to include some additional ad-hoc absorp- 
tion during the calculation of the spectrum (shown as bold 
solid lines in Fig. |2]i. This allows us to produce good agree- 
ment between the synthetic spectra and the Swift observations. 
The resulting total absorption column (including the ISM, the 
CSM from the RG wind and the ad-hoc additional absorption) 



is within the same order of magnitude as the absorpti on col- 
umn required by the fits to the Swift data in iPaper 1 Such 
additional absorption could be due to the presence of a dense 
gas torus clo se to the bin a ry co re, as suggested by the IR ob- 
servations of lEvans et alJ (2007). 

Of course, we have to acknowledge that such a dense torus 
would have an effect on the dynamics of the shocks propa- 
gating through the CSM but this is a problem that can only 
be addressed in 2D or 3D simulations. Furthermore, such a 
dense torus will play a role in the bipolarity seen in the ra- 
dio (O'Brien etal. 2006; Sokoloski et al. 2008; RupenetaH 
2008) and opti cal imagin g ( Bode et al.l 120071 ; iRibeiro et al] 
2002). Indeed, IWalder et al.. (lOOipiave shown in full 3D 
hydrodynamical simulations that the orbiting accreting WD 
leaves a spiral wake in the RG wind along with a slight equa- 
torial density enhancement, and that the high-density accre- 
tion disk around the WD has a substantial size (3x10'^ cm in 
diameter). They illustrate the ejecta col limation efficiency of 
these features during the nova outburst. lOrlando et al.l ( l2009h 
have also computed the X-ray emission from 3D simulations 
of the outburst in a wind with an equatorial density enhance- 
ment but unfortunately do not show the absorption effect of 
the presence of such an enhancement on the resulting syn- 
thetic spectra. 

4.3. Shock position, velocities and flux comparisons 

Figure [3] shows the run 81 forward shock position in milli- 
arcseconds (mas) assuming a distance of 1 .6 kpc as a function 
of time (black solid line), compared to the size measurements 
of three separate features i n resolved radio-inter f erometry im- 
ages. The radio images oflO'Brien et"an (l2008h . iRupen et alj 
(120081) and lSokoloski et al I (l2008h have indeed revealed three 
main synchrotron emitting features (most probably associated 
with shocks) in the remnant of RS Oph: a central expanding 
ring, an eastern lobe further away from the central binary sys- 
tem and an additional western counterpart. The positions of 
the lobes and the size of the expanding ring measured in the 
radio observations are represented by the various points (see 
legend for details). 

It seems natural to compare the position of our spherically 
symmetric shock to the size of the expanding ring. We first 
note that the position of the run 8 1 forward shock is higher 
than the early ring size measurements (t < 30 days) by a fac- 
tor of about two. At later times (30 < f < 50 day s), it is very 
consistent with the ring size measurements of Sok oloski et alj 
(l2008h . before it breaks out of the RG wind (kink circa day 
44). Between days 20 and 45, the model shock position is 
also close to t he po sition of the western lobe measured by 
lO'Brien et alj (120081) . It is however never close to the mea- 
sured positions of the eastern lobe until the shock breakout. 
We conclude that our model shock positions are consistent 
with a distance of 1.6 kpc and that over-estimates at early 
times arise from the fact that we are trying to reproduce the X- 
ray emission of a system with multiple ejection features mov- 
ing at different speeds with a spherically symmetric model. 
The model shock positions sensibly lie between the ring sizes 
and the eastern lobe positions, as a compromise between the 
two. 

Figure |4^ shows the forward shock velocities (solid line) 
measured directly from the hydrodynamical simulation of run 
81 . These are compared to the shock velocities estimated 
in iPaper II from single-temperature fits to the Swift spectra 
(empty squares with error bars) and to the IR velocities from 
iDas. Baneriee & Ashok (2006.) (triangles). The IR velocities 
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Figure 2. Synthetic X-ray spectra for run 81 at days 3, 5, 11, 18, 29, 50, 75 and 100. The thin solid and dashed lines represent a fully neutral and fully ionised 
RG wind, respectively. The bold solid line represents the spectrum including additional absorption. The additional hydrogen equivalent absorption columns are 
(in units of 10^^' cm^"): 4.0 (day 3), 2.5 (day 5), 0.5 (day 11) and 0.3 (day 18). The black crosses represent t he Sw ift data. The dot-dash line in the top left panel 
represents the contribution to the total X-ray spectra from the post-forward shock material (see end of section l43t . 



were taken from the full width zero intensity (FWZl) of the Ol 
and Pa ^ emission lines . The v elocities from the RXTE X-ray 
data of 'Sokolos ki et al.l (l2006h are also shown (crosses). 

Figure ]^ shows the total 1 — 10 keV X-ray fluxes from 
the Swift observations (empty squares) compared to the fluxes 
calculated from run 81 (solid line; including the additional 
ad-hoc absorption). The model fluxes are consistent with the 



observed fluxes. The early turnover after an initial rise around 
day 6 occurs at the same time in both model and observations. 
There is also a clear break in the model curve just after day 
40 which appears to fit the data; this corresponds to the time 
when the model forward shock breaks out of the RG wind. 
This occurs at the same time as the dramatic increase in the 
model velocities. 
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Western lobe (O'Brien et al. 2008; VLBI) 
North-South ring size (O'Brien et al. 2008; VLBI) 




Time (days) 

Figure 3. Run 8 1 forwai'd shock position as a function of t i me (so l id line) compared to the radio imaging measure ments of the size of the expanding ring and the 
position of the eastern and western lobes bv lO'Brien et al.l j2008l) . lRuDen et si\ 120081) and lSokoloski et al] 420081) . The error bars have been omitted for clarity; 
on the VLBI and VLBA measurements they are ~ 1 mas, on the MERLIN data they are ~ 10 mas. 



While the X-ray fluxes are well reproduced, the model 
velocities do not agree with the Paper I velocities derived 
from the X-ray observations. The model velocities are much 
higher, the early turn-over occurs about 4 — 5 days later and 
the deceleration rate is slower. There is however reason- 
able agreeme nt between the v elocit ies obtained in iPaper II 
and those by ISokoloski et al.l (l2006h from RXTE observa- 
tions. It has been pointed out by Tatischeff & Hernanz' (120081) 
and references therein that equation (3) in Paper I for strong 
shock dynamics can underestimate shock velocities where 
non-linear diffusive shock acceleration of particles is efficient, 
which is most probably the case for RS Oph. We also note 
that the immediately post-shock gas is the hottest and would 
provide the best estimate of the shock velocity. However, 
the observed spectrum is made up from contributions from 
a wide range of gas temperatures, emission from the shocked 
ejecta and even perhaps the WD remnant. These will lead to 
a softening of the spectrum and so velocities computed from 
single-temperature fits will be biased towards lower tempera- 
tures and hence lower derived shock velocities. 

The IR velocities shown in Fig. are in better agreement 
with the model forward shock velocities. However it should 
be emphasised that the different methods are likely to be bi- 
ased towards emission from different parts of the expanding 
remnant. This spans a wide range of velocities, including 
the unshocked and shocked fast wind from the WD, and the 
shocked (and unshocked) RG wind. We also show in Fig.|4^ 
the gas velocity just after the reverse shock (blue solid line), 
which is in very good agreement with the IR velocities at early 
times. Taking the FWZI of the IR lines probes only the fastest 
moving material in the system, which suggests that the IR line 



widths are first dominated by the post reverse shock gas and 
only later by the post forward shock gas. 

As previously discussed, in this spherically symmetric 
model the shock wave breaks out of the RG wind just after 
day 40, leading to a dramatic increase in model forward shock 
velocities. By contrast, around the same time (just after day 
30), a decrease in the IR Ol velocities is observed. However, 
our models show that the accelerating breakout material is at 
very low density and also very hot so it will most probably 
not be visible in the IR. After shock breakout, the IR emission 
comes from other regions of the shell where the velocity is 
much lower, thus explaining the strong decrease in IR veloci- 
ties. 

In order to investigate the differences between the directly 
measured shock velocities in our simulations and the veloc- 
ities estimated from the Swift X-ray spectra, we fitted two- 
temperature MEKAL models to our synthetic spectra for run 
8 1 . Single-temperature fits were initially attempted but were 
badly constrained at the high energy end and inconsistent with 
the synthetic spectra. We then estimated forward shock ve- 
locities using the higher of the two temperatures and these are 
also plotted in Fig.|?^ (green solid line). These new velocities 
are in v ery good agreement with the with the velocities de- 
rived bv lPaper H from Swift observations at early times (f < 5 
days) and are much lower than the shock velocities found in 
the hydrodynamic simulation. They then evolve in a similar 
way to the hydrodynamic velocities until around day 20, when 
they start to decrease more slowly than the simulation and the 
velocities derived from Swift observations. Note that the hard 
X-ray count rates in our synthetic model between days 18 and 
50 are slightly greater than the Swift count rates (see Fig.|2|; 
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Figure 4. Run 81 results compared to Swift data from lPaper I and IR data from lDas et aljf2003) . (a) Forward shock velocities: the black empty squares represent 
the Swft data, the red triangles the IR data (filled = Ol; empty = Pa/3) and the solid line the model. The IR velocities were taken from the FWZI/2 of the emission 
lines. The dashed line is the velocities derived from the T fits to our synthetic spectra. The dotted line represents the gas velocity behind the reverse shock. The 
crosses are the Sokoloski et al. 1 2006) RXTE data, (b) X-ray flux (in the range 1 — 10 keV): the empty squares represent the Swft observations and the solid line 
the model. Error bars here are smaller than the size of the squares (< 5%). 



this yields higher temperatures in the fits and consequently 
higher velocities. This explains why the slope of the veloc- 
ity curve derived from our new spectral fits is lower than the 
slope observed in the Swift data during that period. At the 
shock break-out, the computed velocities from the synthetic 
spectral fits show a sharp decrease, much like the IR veloci- 
ties, in contrast to the hydrodynamic velocities which increase 
dramatically, which can be explained in the same way as for 
the IR velocities (see above). 

Figure |5h shows the emission measure as a function of ra- 
dius inside the expanding remnant for run 8 1 at day 3 after 
outburst. As mentioned in section [3] to calculate our syn- 
thetic X-ray spectra we have run through the radial cells of 
our grid binning the emission measures into 30 temperature 
bins evenly (logarithmically) spaced between 10"* — 10^ K. 
The grey scale in Fig. represents the percentage of the total 
emission measure contained within each temperature bin. The 



lower panels show temperature (b), density (c) and velocity 
(d) of the gas. The forward and reverse shocks are clearly vis- 
ible at 6.22 X lO'^ cm and 4.83 x lO'^ cm, respectively. The 
contact discontinuity between the shocked ejecta and shocked 
RG wind (high density, well cooled and therefore low temper- 
ature) is located at 5.8 x lO'-' cm. 

The major contribution to the total emission measure is 
from the region around the contact discontinuity, while the 
broader forward shocked region with temperatures of ^ 10^ 
K (spanning radii fro m 5. 9 — 6.2 x lO'^ cm) contributes only a 
few percent. Orlando et alj (l2009t) also find in their 3D simu- 
lations that most of the X-ray emission originates from a small 
dense region behind the forward shock. This implies that the 
total X-ray emission is dominated by material at lower tem- 
peratures (around the contact discontinuity) and that the con- 
tribution from the hotter forward shocked gas is small. How- 
ever, the shocked gas might still contribute strongly to the 
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Figure 5. (a) Emission measure as a function of radius inside tlie hot sliell for run 8 1 at day 3 after outburst. Tlie data are greyscale coded to represent what 
percentage of EMto( each temperature bin represents. The scaling is shown on the right (ranging 0— 15%). The inset shows the detail of the region around the 
contact discontinuity, which provides the majority of EMtoi. The temperature bin in which the post-shock material lies contributes only a very small fraction to 
EMiot (less than 5%). (b) Gas temperature as a function of radius. The data are greyscale coded to represent which of the 30 temperature bins each part of the hot 
shell belongs to. The lower panels show the gas density (c) and velocity (d) as a function of radius. 



high energy X-rays. 

The dot-dash line in the top left panel of Fig. |2] shows the 
X-ray spectrum emitted by the gas located between the con- 
tact discontinuity and the forward shock (we shall call this the 



post-forward shock gas) with no circumstellar absorption in- 
cluded; the contribution to the total X-ray spectrum is small 
at low energies but is of the order of 50% for energies above 
5 keV. Such a contribution is non-negligible and should have 
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a large impact on the temperature fits to the data. However, 
the temperature of the emitting gas determines the shape of 
an X-ray spectrum while the emission measure and distance 
to the source govern the normalisation of the spectrum. 

The post-forward shock gas spectrum is much flatter than 
the total X-ray spectrum (solid line); by fitting a single tem- 
perature MEKAL model to the post-shock gas spectrum at 
day 3 we obtain ksT = 13.5 keV, translating as a shock ve- 
locity of 3400 km s^^ which is entirely consistent with the 
real shock velocity measured in the simulation (black dots in 
Fig. lit). This flat post-forward shock spectrum gets drowned 
by the bright but cooler emission from the rest of the shell 
which changes the shape of the X-ray spectrum and therefore 
changes the fit temperature. 

Fitting temperatures to the total observed spectra in order 
to estimate post-shock temperatures and therefore shock ve- 
locities is thus unreliable. It also suggests that cosmic ray ac- 
celeration at the shock may not be required to explain the dis- 
crepancies between velocities derived from observed IR line 
widths and those derived from X-ray temperature fits (even 
though such acceleration might still take place). 

Finally, the velocity plot in Fig.|5}l shows that the material 
moving at the fast ejection speed of lO^kms^' contributes 
nothing to the X-ray emission and will thus not be observed, 
while the post reverse shock gas does emit a small amount and 
its velocities are similar to the ones measured by the IR line 
widths (as shown in Fig.H^). 

4.4. The RXTE dataset 

It has to be noted that even though this study primar- 
ily addresses the Swift X-ray data, one cannot ignore the 
higher energy X-ray data obtained with the RXTE observa- 
tory ( Sokoloski et al.l 1200 6). There are indeed some differ- 
ences between the two datasets, taken during the same epoch 
of the remnant's evolution. The break in the Swift data around 
day 7-8 w hich we interpreted as the end of our fast wind phase 
dPaper III) is absent from the RXTE data. Instead, a nearly 
perfect Sedov-Taylor shock deceleration is observed between 
days 3 and 10 (see Fig. 2b in Sokoloski et al. 2006; note that 
once again these velocities were computed from single tem- 
perature fits to X-ray spectra). The discrepancy between the 
data is not fully understood. The RXTE probes a higher en- 
ergy range (3 - 20 keV) than the XRT onboard Swift (0.3 - 10 
keV) which would imply that it is sensitive to only the hottest 
plasma and would therefore give a better estimate of the post- 
shock temperatures and hence the shock velocities. However, 
the observed velocities a re lower than the velocities derived 
from Swift observations ( IPaper J) which is inconsistent with 
our previous analysis showing that the single temperature fits 
to the Swift data underestimate the postshock temperatures in 
the system. 

Another important difference between the datasets is that 
while an increase in total integrated flux is reported for the 
Swift data between days 3 and 5, a decrease is observed by 
the RXTE. This can be explained by the fact that a large por- 
tion of the energy radiated at day 3 is in the high energy range 
(thus not well detected by Swift) and then as the blastwave 
cools adiabatically the radiation moves to lower energies re- 
sulting at day 5 in an increase of 5 — 10 keV photon counts for 
Swift and a decrease of 7 — 20 keV photon counts for RXTE. 
In addition, as the forward shock progresses outward clearing 
the RG wind, the amount of circumstellar absorption reduces 
which also contributes to an increase of soft X-ray flux. An 



increase of 3 — 5 keV photon counts is observed in the RXTE 
data which confirms this statement. However, if one compares 
the same energy band in the two datasets between days 3 and 
5, one can see that the X-ray counts between 5 — 10 keV in- 
creases in the Swift data while it decreases in the RXTE data. 
Inconsistent datasets will inevitably be difficult to fit with a 
single model. 

The reason why a Sedov-Taylor model fits the early RXTE 
data and not the Swift is unclear. A general Sedov-Taylor 
blastwave is produced by an instantaneous release of energy 
at a point inside an ambient medium with a power-law den- 
sity profile. No mass is ejected, there is no reverse shock, 
and the expanding shell is made up only from the ambient 
gas swept up by the shock. In our model, we have adopted 
a more physically realistic scenario likely to be applicable to 
the case of RS Oph where some mass has been ejected in the 
outburst. Here, the outburst energy is delivered via the kinetic 
energy of a wind which then interacts with a circumstellar 
medium with a density described by an inverse-square law 
as produced by a stellar wind. We have experimented with 
varying the amounts of outburst energy in kinetic energy or 
in thermal energy. For example, we have run a simulation 
using the parameters from run 2 (same Mgj and Eq) but with 
half of the outburst energy in thermal form and half in kinetic 
form. The results are quite different at early times, while re- 
maining very similar at later times (f > 30 days). At early 
times, the blastwave expands much faster (due to the instan- 
taneous release of half the outburst energy), the X-ray spectra 
are much flatter (yielding higher X-ray temperatures) and the 
velocities are much more Sedov-like. As opposed to the ini- 
tial run 2, there is a decrease of X-ray flux between days 3 and 
5 when folded through both the Swift instrument response and 
the RXTE response. This is more in line with the RXTE ob- 
servations but no longer in agreement with the Swift data. A 
more extreme example injecting 100% o f the outburst en ergy 
in thermal form (as in the models of O' Brien et alj|1992l) also 
shows a decrease in 5 — 10 keV photon counts when folded 
through the Swift response. We have not been able to obtain 
results which are consistent with both the Swift and RXTE 
observations. 

5. CONCLUSIONS 

We have t aken the r esults from our hydrodynamical models 
described in 'Paper II and calculated synthetic X-ra y spectra 
using the XSPEC. The results from runs 1 to 16 of iPaper 13 
were first compared to the Swift spectra. It was found that the 
rate of decrease of the hard X-ray counts was related to the 
rate of decrease of the forward shock velocities, and a best 
fitting model would need a relatively low ejecta mass. How- 
ever, runs with very high ejection velocities V2 overestimated 
the X-ray counts across the spectrum at all times. Run 13 
which was thought to be the best fit (based on the forward 
shock velocities) in Paper III was thus excluded. 

A wider exploration of the ID model parameter space was 
undertaken to find a model which best fit the observed spec- 
tra. A total of 88 runs were carried out, and fluxes were ex- 
tracted for each day between day 1 and day 100 after out- 
burst for each run. These model fluxes were compared to the 
fluxes from 70 Swift observational epochs during the same 
period, being careful to exclude SSS emission as appropri- 
ate. Large differences between the observed and theoretical 
fluxes allowed us to exclude many runs, and the spectra of the 
few remaining runs with good matching fluxes were exam- 
ined in order to find a best-fitting model. We are now confi- 
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dent that the mass ejected in the RS Oph outburst is relatively 
low, around (2 — 5) x 10^^ M©, and that the shock veloc- 
ities are higher than suggested by X-ray single temperature 
fits. It was however never possible to simultaneously obtain 
a sufficiently high level of hard X-ray emission, a flat mid- 
energy spectrum and an appropriately low level of emission 
at low energies without including some additional ad-hoc ab- 
sorption, thus making the derived densities in the RG wind 
uncertain. 

The fast wind velocity of the best-fitting model 
(10000 km s^') is considered unusu ally high for a nova 
out burst, but radio obser vations of lO'Brien et al.l (l2006h 
and ISokoloski et al.l (l2008h show evidence for bipolar lobes 
moving at speeds of the order of 5000 kr n s~' at earl y times 
(f < 50 days). Chandra observations of ILuna et al.i (l2009l) 
also repor t a jet moving at sp eeds greater than 6000 kms^^ 
Similarly. iRibeiro et al.l (120091) deduce high ejection veloci- 
ties of - 5100 kms"' from HST imaging and ground-based 
optical spectra, albeit the latter two studies are based on 
observations taken much later (538 and 155 days after 
outburst, respectively) than the timescales addressed in this 
paper. Furthermore, the fast ejection velocity may not be 
easily directly observed as the ejecta are decelerated when 
they run into the circumstellar medium. 

Another important conclusion we can draw from this study 
is that in the Swift/XRT energy band, the emission from the 
post-shock material is negligible compared to that of the re- 
gions around the high-density contact discontinuity, and that 
fitting emission models with single temperatures (even when 
several components with different temperatures are included) 
to X-ray spectra greatly underestimates immediate post-shock 
temperatures and hence forward shock velocities. We also 
showed that velocity estimates from line widths (for example 
the IR line widths in the case of RS Oph) also underestimate 
the forward shock velocities. 

Some inconsistencies are however present between the 
Swift and RXTE X-ray datasets, in particular in the forward 
shock velocities between days 3 and 10. While the spectra of 
our best fitting model agree well with the Swift observations, 
they do not fit the RXTE data. This does not affect our results 
for the subsequent evolution of the remnant (in particular the 
ejected mass) since RS Oph was very faint in the RXTE en- 
ergy range after day 10, as opposed to Swift. It just means that 
the manner in which the mass and energy are ejected during 
the outburst might be different (for example, the proportions 
of outburst energy in the form of kinetic and thermal energy 
may be different). 

It has proved difficult to find a set of parameter values for 
this spherically symmetric model which matches the observed 
X-ray spectra for all epochs during the first 100 days of the 
remnant's evolution. However, radio and optical observations 
have shown that the remnant of RS Oph has a bipolar mor- 
phology. Our suggestion that the additional absorption re- 
quired to fit the low-energy X-ray spectra might be present 
in the form of a torus would also fit with an evolving bipolar 
morphology. Preliminary work on 2D codes with a bipolar 
ejection is underway, using an equatorial density enhance- 
ment in the RG wind. Initial results have shown that is it 
possible to obtain fast moving material in the polar direction 
yielding hard X-rays and a high absorption column along the 
equator at the same time, better res embling the observed X- 
ray spectra of RS Oph ( Vavtetll2009 ). A bipolar model might 
also be able to fit the RXTE observations better if the remnant 



exhibits more of a Sedov-Taylor blastwave behaviour along 
one direction. Much work however remains, due to the in- 
creasingly large parameter space worthy of exploration. 

The authors are very grateful to the Swift Mission Op- 
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edges support from research grant ANR-06-C1S6-009-01 for 
the SiNeRGHy project at the CEA/SAp Saclay. The authors 
would also like to thank the referee J. Sokoloski for very use- 
ful comments which have helped to improve this paper and 
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radio observations in Fig. [3] 
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